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Abstract 

The  field  diffracted  by  an  aperture  of  any  shape  is  evaluated  asymptotically 
for  small  wavelengths  on  the  basis  of  several  approximate  diffraction  theories. 
These  are  the  Kirchhoff  method,  the  two  customaiy  modifications  of  it,  and  W, 
Brauhbek's  new  modification  of  this  method.  In  each  case  a  double  integral  over 
the  apertiu^  is  evaluated  asjTnptotically,  and  contributions  from  interior  stationary 
points,  edge  stationary  points,  and  comers  of  the  edge  are  obtained.  The  contribu- 
tions from  points  of  each  type,  according  to  all  the  theories  examined,  are  exactly 
of  the  form  recently  deduced  by  J.  B.  Keller  using  his  geometrical  theory  of  diffrac- 
tion. The  edge  stationary  points  and  the  comers  correspond  respectively  to  rays 
singly  diffracted  from  edges  and  from  comers.  For  edge  diffracted  rays  the 
Braunbek  diffraction  coefficients  alone  coincide  with  those  given  by  the  geometrical 
theory  of  diffraction,  but  the  Braunbek  method  does  not  apply  to  comer  diffraction. 
None  of  these  methods  takes  account  of  multiple  diffraction  as  does  the  geometrical 
theory  of  diffraction. 
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1.  Introduction 

Diffraction  by  an  aperture  in  a  thin  screen  was  treated  by  a  new  method 
in  Part  I"-  ■'  of  this  paper.  In  this  second  part  we  compare  this  method  with 
certain  other  methods  which  have  been  used  in  the  past  to  treat  the  same  problem. 
The  other  methods  which  we  consider  are  the  Kirchhoff  method,  the  two  customary 
modifications  of  this  method'-  -'  and  W.  Braunbek's  new  modification  of  it*-  •'• 
These  methods  have  been  carefully  discussed  by  C.  J,  Bou^rfcamp"-  -* , 

The  geometrical  theory  of  diffraction  used  in  Part  I  involves  the  concept 
of  a  diffracted  ray.  It  enables  us  to  express  the  singly  diffracted  field 
explicitly  in  terms  of  geometrical  quantities.  Similar  geometrical  expressions 
can  also  be  obtained  for  the  multiply  diffracted  field.  In  the  Kirclihoff  method 
and  its  modifications,  on  the  other  hand,  the  diffracted  field  is  expressed  as  a 
double  integral  over  the  aperture,  V7e  vrill  evaluate  the  double  integrals  given 
by  all  these  methods  asymptotically  for  small  wavelengths  X.  We  will  find  that 
eacii  one  leads  to  the  formula  deduced  in  Part  I  for  the  singly  edge-diffracted  field, 
with  an  appropriate  diffraction  coefficient  in  each  case.  This  fact  supports  the  theory 
used  in  Part  I.  Only  the  Braunbek  method  yields  exactly  the  same  edge  diffraction 
coefficient,  and  therefore  the  same  singly  edge-diffracted  field,  as  is  obtained  in 
Part  I,  This  agreement  between  the  Braunbek  method  and  the  geometrical  theory  of 
diffraction  lends  support  to  both  of  these  methods,  and  shows  that  the  simpler 
geometrical  theory  can  be  used  instead  of  the  more  complicated  Braunbek  method. 

Neither  the  Braunbek  method  nor  the  Kirchhoff  methods  take  account  of 
multiple  diffraction,  which  is  included  in  the  theory  of  Part  I.  Furthermore,  the 
Braunbek  method  cannot  be  applied  to  apertures  with  comers.  The  other  Kirchhoff 
methods  can  be  applied  to  such  apertures  and  yield  results  of  the  form  given  in 
Part  I  for  the  field  on  a  ray  diffracted  from  a  comer.  However  the  Kirchhoff 
diffraction  coefficients  differ  from  those  of  Part  I,  which  come  from  exact 
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solutions  of  certain  diffraction  problems  involving  comers.  Therefore  the 
Kirchhoff  diffraction  coefficients  are  inaccurate,  although  they  are  easy  to 
compute.  It  is  possible  to  modify  the  Braunbek  method  so  that  it  applies  to 
apertures  with  comers  and  yields  results  for  them  in  agreement  with  those  of 
Part  I.  But  even  then  it  will  not  take  account  of  multiple  diffraction. 

When  the  transmission  cross  section  of  an  aperture  is  computed  by  any  of 
the  Kirchhoff  methods  using  the  cross  section  theorem,  only  the  geometrical  optics 
cross  section  is  obtained.  The  Braunbek  method  yields  a  correction  of  oixier  k" 
to  the  geometrical  cross  section.  However  the  geometrical  theory  of  diffraction 
yields  a  correction  of  order  k"^'^  for  a  circular  hole  in  a  screen  on  which  u<.  The 
failure  of  the  Kirchhoff  and  Braimbek  methods  to  yield  the  proper  correction  in  this 
case  can  be  attributed  to  their  failure  to  account  for  multiple  diffraction. 

In  spite  of  the  defects  of  the  Braunbek  and  Kirchhoff  methods  pointed  out 
above,  they  do  have  one  advantage  over  the  method  of  Part  I.  That  method  yields 
infinite  values  on  caustics  of  the  diffracted  rays  and  on  shadow  boundaries.  It 
must  be  supplemented  by  special  considerations  in  order  to  yield  finite  results 
at  such  places.  The  Braimbek  and  Kirchhoff  methods  do  not  have  this  defect. 
Therefor©  they  can  be  used  to  supplement  the  method  of  Part  I  on  the  surfaces 
where  it  fails. 

For  the  special  case  of  an  incident  spherical  wave.  A,  Rubinowicz"-  J  has 
evaluated  the  original  Kirchhoff  integral  asymptotically.  He  did  it  by  transforming 
the  double  integral  into  a  line  integral  along  the  rim  of  the  aperture  and  then 
evaluating  the  line  integral  by  stationary  phase.  His  result,  which  he  interpreted 
in  terms  of  diffracted  rays,  is  a  speci<il  case  of  the  resvdt  which  we  will  obtain. 
The  integrals  given  by  the  modified  forms  of  the  "^irchhoff  method  were  evaluated 
asymptotically  by  N.G.";''an  Kampen'-  ■' .  His  results  are  the  sane  as  those  we  obtain 


in  ohese  cases  and  his  heuristic  method  of  evaluation  has  been  justified  by 
J,  Focke*-  -'  and  N,  Chakol-'J,  However  neither  of  their  analyses  is  directly 
applicable  to  the  integrals  which  occur  in  the  Braunbek  method.  These  integrals 
are  our  main  concern,  since  they  have  not  previously  been  evaluated  asymptot- 
ically. We  have  treated  them  by  a  method  different  from  those  given  in  [5],  [6], 
[7] .  For  the  sake  of  xmiformity  we  have  also  treated  the  modified  Kirchhoff 
integrals  by  this  same  method. 


Figure  1 
The  coordinates  cand  s  of  a  point  in  the 
aperture.  Distance  along  the  edge  from  some 
fixed  point  is  denoted  by  ct;  The  cr  coordinate 
of  a  point  is  the  value  of  o" at  the  nearest 
point  on  the  edge.  The  s  coordinate  is  the 
distance  from  the  point  to  the  edge. 
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2.  The  Kirchhof f  and  Braunbek  methods 

Consider  a  thin  screen  S  containing  an  aperture  H  and  l3rLng  in  the 

plane  z  =  0,  Suppose  an  incident  scalar  wave  u.   comes  from  the  left  (z  <  O) 
^^  mc 

and  strikes  the  screen.  We  wish  to  determine  the  resulting  field  u,  assuming 

that  it  satisfies  the  reduced  wave  equation.   The  total  field  u  must  consist  of 

the  incident  field  u.   ,  to  the  left  of  the  screen,  and  an  additional  scattered 

mc 

or  outgoing  part  on  both  sides  of  the  screen.  Furthermore  on  the  screen  either 
u  or  its  normal  derivative  must  vanish,  depending  upon  the  type  of  screen. 
Thus  u  satisfies  one  of  the  boundary  conditions 


(1)  u  =  0       on  S 

(2)  g  =  0      onS. 


By  means  of  Green's  theorem,  the  field  u  behind  the  screen  (z  >  O)  can 
be  represented  as  an  integral  over  the  plane  z  =  0  in  three  different  ways: 


(3)  u(P)       ^ 


H^ 


r  r       -,    /  ikR\         xkR  ^ 
5    (  e )       e Su 

^  8z  \  R      /    ~       R     az 


djcdy, 


"•^        -"W'^l  -k  ['^)^^. 


*''     "(« =  -  5i;  f  ^  i  *^^- 


In  these  equations  R  denotes  the  distance  from  the  point  P,   at  which  u  is 
evaluated,   to  a  point  on  the  plane  z  =  0.     In  the  integrands,  u  and  ^  are  to 
be  evaluated  on  the  back  side  (z  =0+)  of  the  screen.     The  propagation  constant 
k  =  2TiA,>fhere  \  is  the  wavelength. 

In  order  to  compute  u(P)  from  the  first  of  these  fomailas  we  must  know 
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u  and  —  on  the  plane  z  =  0,  Vie   know  that  one  or  the  other  of  these  quantities 

is  zei^3  on  S  according  as  bovmdary  condition  (l)  or  (2)  is  satisfied,  Kirchhoff 

assumed  that  both  u  and  rr—  are  zero  on  the  rear  of  S  and  that  u  =  u.  %  ■k-  "  — =? 

dz  inc'  oz    dz 

in  the  aperture  H.  Then  (3)  becomes 

T  r  r      t>  /  ikR  \     ikR   Su.  -\ 

(«'      V«  =  i J,  kc  * (V)  -  V  ^ ]  ^. 

We  will  call  the  field  ^en  by  this  method  the  Kirchhoff  approximation  and  denote 
it  by  Ujj(P). 

If  (1|)  is  used  instead  of  (3),  only  the  value  of  u  need  be  known. 
Suppose  that  (l)  is  satisfied.  Then,  if  we  assume  that  u  =  u.   in  H,  (U)  becomes 


inc 


(7)      u,(P)=:^|  u,.„  I- (41^ )  cbcdy. 


V^'--Ji;j   "mo  ^iVj 


H 

We  will  call  the  field  given  by  (?)  the  first  modified  Kirchhoff  approximation 
and  denote  it  by  U-(P).  In  a  similar  way,  if  (5)  is  used  and  (2)  is  satisfied,  the 

au. 

assumption  "5—  =  -^ in  H  leads  to 


T  /■      ikR   8u. 


We  will  call  the  field  given  by  (8)  the  second  modified  Kirchhoff  approximation 
and  denote  it  by  Up(P). 

The  three  fields  Uj^,  u-  and  Up  are  all  approximate  descriptions  of  the 
diffracted  field  since  they  are  based  upon  approximate  values  of  u  and/or  t-  on 
the  plane  z  =  0.  Physically  it  seems  clear  that  the  assumed  values  on  this  plane 
should  be  the  more  nearly  correct  the  smaller  the  wavelength  X  or  the  larger 
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the  value  of  k.     Therefore  we  will  evaluate  the  integrals  which  represent  these 
fields  asymptotically  for  k  large.     Before  doing  so  it  is  convariient  to  notice  that 

(9)  ^  °  7  (^1  *  ^2^* 

Therefore  ire  need  only  evaluate  u,  and  Up,  and  then  we  can  find  u^   from  (5). 

The  evaluations  of  the  integrals  are  contained  in  the  next  sections. 

The  new  modification  of  the  Kirchhoff  method  proposed  by  Braunbek 

makes  use  of  the  representation  (U)  when  u  =  0  on  S,  and  of  (5)  when  q^  =  0  on  S. 

In  the  case  u  =  0  on  S  it  assum.es  that  u  =  u.  „  in  H  except  near  the  edge.  Near 

any  point  Q  on  the  edge  u  is  assumed  to  be  given  by  the  field  near  a  half-plane 

on  which  u  =  0.  The  half- plane  lies  in  the  plane  of  the  screen,  and  its  edge 

is  tangent  to  the  edge  of  the  screen  at  Q.  /  plane  wave  is  incident  on  it  from 

the  direction,  and  with  the  amplitude,  of  the  field  incident  on  the  aperture 

at  Q.  Since  the  field  diffracted  by  a  half-plane  is  known,  an  explicit 

expression  for  the  field  near  any  point  Q  on  the  edge  is  available.  In  the 

case  |i^  =  0  on  S  it  is  assumed  that  -^  =   ^°  in  H  except  near  the  edge.  Near 
6z  az       az 


the  edge  ~  haa  the  yalue  given  by  the  field  near  a  half- plane  on  which  ^r  ■  0. 
We  will  call  the  two  fields  given  by  (U)  and  (5),  with  the  Braunbek  integrands, 
n^  and  tUg  respectively. 

Just  as  in  the  case  of  the  Kirchhoff  approximation  and  its  usual 
modifications,  it  is  to  be  expected  that  the  Braunbek  approximation  becomes 
more  accurate  as  k  increases.  Therefore  we  will  also  evaluate  the  integrals 
given  by  this  approximation  for  large  k.  It  also  seems  likely  that  this  approxi- 
mation is  more  accurate  than  the  previous  ones.  However  it  is  clear  that  the 
Braunbek  method  does  not  apply  near  comers  of  the  edge.  The  obvious  modifica- 
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tion  required  at  such  comers  is  that  the  field  there  should  be  replaced  by 

the  field  near  an  'infinite  comer,'  i.e.  a  screen  with  two  straight  edges  meeting 

at  a  comer. 

In  evaluating  the  Kirchhoff  and  Braunbek  integrals  we  will  assume  that  the 

incident  field  has  the  foim 
(10)      u^^  ..  Ae^ 

The  phase  T  and  amplitude  A  are  functions  of  position.  Since  a  field  of  this 
foni  is  not  necessarily  a  solution  of  the  reduced  wave  equation,  we  have  used 
the  sign  ''~'  in  (10).  This  indicates  that  the  incident  field  is  only  asymptoti- 
cally of  this  form,  for  large  k. 

3.  Asymptotic  evaluation  of  the  Kirchhoff  integrals 

Upon  inserting  the  incident  field  (10)  into  (?)  and  (8)  we  obtain 


(n)     u^(P)^  §  I  AR-\  e^(^*^  W, 

H 

(12)     U2(P)  -  -  i  j  AR-^,  e^^^^^W. 


H 


Note  that  in  both  equations  only  the  leading  terms  in  k  have  been  retained. 

When  an  integral  such  as  that  in  (ll)  or  (12)  is  evaluated  asymptotically 
for  k  large,  the  entire  contribution  to  the  asymptotic  value  comes  from  the 
neighborhoods  of  certain  critical  points L  J»l'-i,  Ihese  points  are  of  three  kinds. 
Critical  points  of  the  first  kind,  which  may  be  called  interior  stationary  points,  sue 
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interior  points  of  II  at  which  the  phase  R+T  is  stationary.     Critical  points 

of  the  second  kind,  which  may  be  called  edge  stationary  points,  are  points  of 

the  edge  at  which  the  phase  is  stationary  vrLth  respect  to  variations  along  the 
edge.     Critical  points  of  the  third  kind,  or  comers,  are,  to  the  order  we  consider, 

points  of  the  edge  at  which  the  slope  is  discontinuous.     Let  us  designate  the 

cont:lbutions  of  the  three  kinds  of  stationary  points  to  u^(P)  by  u.  (P),  u.    (?)  and 

uT     (?)  with  a  similar  notation  for  \X2(?).     Then  we  have 


(13)      U3_(P)  ^  uJ(P)  +  uJ^(P)  +  uJ^^(P), 

(Hi)      U2(P)  ~  uJ(P)  +  u^^(P)  +  uJ^^(P). 

2       2 

At  an  interior  stationary  point  (R+T)  =  (R+l)  =  0.  Since  (VT)  =  (VR)  =  1, 

X       y 

while  T-  >  0  and  ^z  *^  °»  ^*  follows  that  VY  =  -VR.  Thus,  noting  that  VR  points  from  ? 
toward  the  aperture,  we  see  that  the  incident  ray  in  the  direction  VT  must,  when 
extended,  pass  through  P.  Thei-efore  to  locate  critical  points  of  the  first  kind  we 
first  extend  all  the  incident  rays  through  the  plane  z  •  0.  If  a  ray  which  passes 
through  P  goes  through  the  aperture  H,  then  its  point  of  intersection  with  H  is  a 
critical  point  of  the  first  kind.  Thus  if  ?  is  in  the  geometrically  lit  region,  there 
is  such  a  stationary  point,  while  if  P  is  in  the  geometrically  dark  region  there 
is  none. 

The  contributions  u,(P)  and  U2(P)  of  the  neighborhoods  of  interior 
stationary  points  to  the  integrals  in  (ll)  and  (12)  are,   up   to  terms  of 
order  k"  , 

T      T         .-4  ▼     ik(R+T)     , 
(15)      u^^P)  =  uJ(P)  -  21      ^      ^ +  0(iw). 
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Here  e  and  A  are  defined  by 


(16)     A  =  (R*^W^*^^yy  -  ^^"-^V 


^+1   if  A  >  0  and  (R+Y)_  >  0 


(17)      e=< 


XX 

-1   if  A  >  0  and  (R+Y)^  <  0 
-i   if  A  <  0  , 


v» 


All  quantities  in  each  term  in  (l5)  are  to  be  evaluated  at  an  interior 

stationary  point,  and  the  sum  extends  over  all  of  these  points.  Since  at  such 

points T  "-R., both  (ll)  and  (12)  yield  the  same  result.  Equation  {!$)   can  be 
z   z 

obtained  from  (ll)  or  (12)  by  expanding  the  phase  of  the  integrand  in  a  Taylor 
series  (up  to  quadratic  terms)  around  each  stationary  point  and  then  appljring 
the  usual  method  of  stationary  phase  twice.  It  can  also  be  obtained  by  making 
use  of  equation  (109  )  of  J,  Focke'-  ■'.  If  ^«  O  at  any  point,  then  P  lies  on  a 
caustic  of  the  incident  rays  and  that  point  contributes  a  more  complicated  term  to  (13 )< 
In  Appendix  I  it  is  shown  that  (l$.)  is  exactly  equal  to  u.  -(P). 

Therefore  we  define  u    (P)  as  zero  if  P  is  in  the  shadow  and  as  u.   (P)  i^  P 

geom  mc 

is  in  the  lit  region.  Then  we  may  say  that 


(18)      ^(P)  =  uI(P)  .  u   (P). 


If  P  lies  on  a  shadow  boundaiy  the  corresponding  interior  stationary  point  becomes 
a  boundary  stationary  point.     In  this  case  only  half  the  neighborhood  lies  in  the 
aperture,  provided  that  the  stationary  point  is  not  at  a  comer.     Thus  the  contid- 
bution  of  this  point  is  just  one-half  of  the  typical  term  in  (l5).     If  there  is 
only  one  stationary  point  then  that  point  contributes  i  u.     (P). 

b        J.T1C 
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Let  us  now  consider  boundary  stationary  points.      We  could  obtain  the 
contributions  from  such  points  by  using  Focke's  equation  (  123    ).     However  it 
is  simpler  to  evaluate  these  contributions  directly.     To  this  end  we  introduce 
the  arclength  a    along  the  edge  and  the  distance  s  along  the  normal  to  the  edge  as 
new  variables  (see  Figure  1),     In  a  sufficiently  narrow  strip  along  any  smooth  portion 
of  the  edge  each  point  has  just  one  pair  of  coordinate  vaDues  7-    and  s.     The  area 
element  in  these  coordinates  is    \l  -  sp~  (cr)|do-ds  where  p(o-  )   is  the  radius  of 
curvature  of  the  edge,  measured  positively  when  the  center  of  curvature  is  on  the 
same  side  of  the  edge  as  the  aperture.     In  terms  of  these  coordinates  we  have 


(19)  u^^(P)    -  g  [  (pjC\  [l  .  p-^(a-)s]e^(^*^^dsdcr-. 

In  (19)  L  denotes  the  length  of  the  edge.  The  tipper  limit  on  the  s  integration 
has  not  been  indicated  since  the  only  contributions  to  u.  (P)  come  from  the 
edge,  which  is  at  s  •=  0, 

We  now  integrate  by  parts  with  respect  to  s  in  (19)  and  retain  only 
the  leading  term  in  k.  We  also  ignore  the  contribution  from  the  upper  limit 
in  accordance  with  the  definition  of  a,'^(P),  Then  (19)  becomes 

The  integral  in  (20)  can  be  evaluated  by  stationary  phase  and  yields 


(21)      ^^h?)   ^  Z  V^'^^^*''^^  • 
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The  quantities  G  and  IL   in  (21)  are  defined  by 
(22)  G  =  — §2£J e'^'^^  -  ^"^, 


_lnA 

(23)         n 


-e  Rj5_ 


^        /2iS  sin  p(T  +R  )  ' 

e  9       0 


In  (21)  the  sum  extends  over  the  edge  stationary  points,  at  each  of  which 
T_  +  R  =0.  Each  terra  in  the  sum  is  evaluated  at  one  of  these  points.  The 

O"    O" 

angle  p  at  each  point  is  the  angle  between  VT  (the  incident  ray  direction)  and 
the  edge.  The  sign  in  the  eocponent  of  (22)  is,  at  each  point,  the  same  as  the 
sign  of  (T+R)^^. 

The  condition  for  a  point  on  the  edge  to  be  an  edge  stationary  point 
is  T^  +  R^=  0.  This  can  be  written  as  Vr»t  "  -  VR*t,  where  t  is  a  unit  tangent 
to  the  edge.  In  the  latter  form  this  condition  can  be  interpreted  geometrically. 
It  states  that  the  incident  ray  makes  the  same  angle  with  the  edge  as  does  the 
line  from  the  edge  to  P.  Furtbermoi^  the  incident  ray  and  the  line  to  P  lie 
on  opposite  sides  of  the  plane  normal  to  the  edge  at  the  stationary  point. 
If  we  call  the  line  from  the  edge  to  P  a  diffracted  ray,  this  condition  is 
just  the  law  of  diffraction''  -'.  Thtis  we  see  that  tu  (?)  contains  exactly  one 
term  for  each  edge  diffracted  ray  through  P,  and  no  other  terms. 

Precisely  the  same  evaluation  applies  to  Up  (P)  if  we  replace  -R  by 
T  in  the  integrand.  Thus  the  final  result  for  \u   (P)  is  also  given  by  (21), 
but  with  D^  replaced  by 

(21;)     Dp 2 la ^ 

'^       yZiSi  sin  p(T  +R  ) 

3   S 
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Then  from  (9),  tij,  (P)  is  also  given  by  (21 )  with  D-  replaced  by 

(25)  Dj^  =  I  (Dj^+Dg). 

In  Appendix  II  it  is  shown  that 

(26)  G-  [R(l+f^^)]-^/2. 

Here  R  denotes  the  distance  to  P  along  the  diffracted  ray  from  the  point  of 
diffraction  and  -p,  denotes  the  distance  to  the  caustic  of  the  diffracted  rays. 
With  this  value  of  G  inserted,  (21)  is  exactly  the  same  as  equation  (l?)  of 
Part  I.  That  equation  expresses  the  field  at  P  on  all  rays  singly  diffracted 
from  the  aperture  edge.  Therefore  ru  (P)  may  be  interpreted  as  the  contribution 
to  u- (P)  from  such  rays*  Then  D,  is  the  diffraction  coefficient  for  edge 
diffraction  given  by  the  first  modified  Kirchhoff  approximation.  Since  u-  (P) 
and  u~  (P)  are  also  given  by  (21)  they  also  correspond  to  the  field  on  rays 
singly  diffracted  from  the  edge,  D  and  Du  are  then  the  edge  diffraction 
coefficients  given  by  the  second  modified  Kirchhoff  approximation  and  the  usual 
Kirchhoff  approximation  respectively. 

The  diffraction  coefficients  D- ,  Dp,  and  Pj,  can  be  written  in  a  more 
geometrical  form  in  terms  of  the  angles  a,  9  and  p  introduced  in  Part  I.  As 
va  have  already  explained,  p  is  the  angle  between  Vf  (the  incident  ray)  and  the 
edge.  To  define  a  and  9  we  project  the  incident  and  diffracted  rays  (VI  and  -VR) 
into  the  plane  normal  to  the  edge  at  the  edge  stationary  point  (the  point  of 
diffraction).  The  angles  between  these  projections  and  the  neg:ative  z-dii«ction 
are  a  and  0  respectively  (see  Figure  2),  In  terms  of  these  angles  we  have 
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Incident  Roy 


Diffracted 
Ray 


Screen 


Figure  2 
The  projections  of  the  incident 
and  diffracted  rays  into  a  plane 
normal  to  the  edge  at  the  stationary- 
point.  The  angles  a  and  &  are  the 
angles  between  these  projections 
and  the  negative  z  axis,  measured 
as  sbovm  in  the  figure.  The  edge 
is  perpendiciJ.ar  to  the  plane  of 
the  figure. 


Figure  3 
The  X,  y  coordinates  at  a  corner 
of  the  edge.  The  angle  between 
the  tangents  to  the  edge  at  the 
corner  is  d. 
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R  «  sin  p  cos  9,  T  ■  sin  ,8  cos  o,  R  ■  -sin  p  sin  6  and  T  «=  -sin  6  sin  a. 
z  z  s  s 

Inserting  these  expressions  into  (23)>  (2U),  and  (25)  yields 
(27)     Ilj^  -    *    °°^  ® 


/2jTk  sin  p(sin  9  +  sin  a) 


inA 
(28)  D,    -     -^  °°^^ 


/?»?<  sin  p(sin  9  +  sin  a) 


*  A 

(29)     D^  =   e     (cos  9  -  cos  a) 

2>^Tnc  sin  p(sin  9  +  sin  a) 

It  is  to  be  noted  that  these  diffraction  coefficients  become  infinite 

trtien  9  =  n  +  a,  which  occurs  if  P  lies  on  the  shadow  boundary.  This  is  due  to 

the  vanishing  of  the  factor  Y  +  R  in  the  denominators  of  (23),(2ii)  and  (25). 

s        s 

When  this  factor  vanishes,  the  phase  at  the  edge  stationary  point  is  also 
stationary  with  respect  to  variations  into  the  interior.     The  integration  by 
parts  which  led  to  (20)  is  then  invalid  and  another  method  of  evaluating  the 
integral  (19)  is  required.     The  resxilt  of  such  an  evaluation  was  referred  to 
after  equation  (18)  anj  will  not  be  considered  further. 

Finally  we  must  consider  u,      (P),   the  field  contribution  from  the  comers. 
At  each  of  these  points  the  tangent  to  the  edge  is  discontinuous,  and  we  will 
denote  by  d  the  angle  between  the  tangents  at  a  typical  comer  (see  Figure  3). 
Then  if  the  x-axis  is  tangent  to  one  edge  at  this  comer  and  if  the  origin  is  at 
the  vertex,  the  other  edge  is  tangent  to  the  line  y*  x  tan  d.     We  will  write  the 
equations  of  the  two  edges  as  y  =  y,(x)  and  y  ■=  y2(3c),  where  yJ!(0)  ■  0  and 
y'(0)   =  tan  d.     Now  we  can  write  u,      (P)   in  the  form 
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(30) 


uJ^^(P) 


AR-V,"'(''*»)dydx. 


The  sum  extends  over  the  comers  of  the  edge  and  at  each  corner  the  coordinates 
are  situated  in  the  manner  just  described.  The  upper  limit  of  the  x  integration 
has  been  omitted,  since  a,  (P)  includes  only  the  contribution  to  the  integral 
from  the  comer,  which  is  at  the  lower  limit  x  »  0, 

We  now  integrate  by  parts  in  (30),  first  with  respect  to  y  and  then  with 
respect  to  x.  Each  time  we  retain  only  the  leading  term  in  k,  which  comes  from 
the  endpoints.  In  this  way  we  obtain 


(3X,    „J"(P,.Z||f«-\|^ 


dx 


jy^Cx) 


^  2nik(R+T; 


[R(x,y^)+T(x,y^)]^"  [R(x,y2)+Y(x,y2)]. 


This  can  be  rewritten  in  the  form 


(32)      u^^(p)  ^  r  V  ^ 


ik(Y+R) 


wheiTe 


(33) 


R 
1         z 


"5nk    Tr+t7 


1  1 , 

|.RU,yi  ^*^^^»yi^jx  "     p(x,y2)+T(x,y2)J^ 


-  Hi  - 


In  order  to  express  c^   in  geometrical  terms,  we  introducse  the  angles  a  and  b 

betvreen  VI  and  the  two  edges  at  the  comer,   and  the  angle  c  between  VT  and  the 

z-axLs  (see  Figure  U).     Similarly  we  define  a'  and  b' as  the  angles  between  -VR  and  the 

two  edges,   and  c'   as  the  angle  between  -7R  and  the  z-axis.     Then  R     =  -  cos  c', 

R^+I^  ■  cos  a  -  cos  a«,and   [R(x,y2)+T(x,y-)j     -  cos  b  -  cos  b'.     Upon  using  these 

cos  d 
relations  in  (33)  we  obtain 


/-I  \    c  ■  .\       cos  c*  sin  d 

1  '  2nk  (cos  a-cos  a'Kcos  b-coy  b'^ 


A  similar  evaluation  shows  that  u^  (P)  is  also  given  by  (32),  with  C^ 
replaced  by  C-.  We  obtain  C«  by  replacing  -R  by  T  in  (33),  Thus  we  find 

/-^\    r  ■  \       cos  c  sin  d 

■^  '     2  "  2i3c  (cos  a-cos  a') (cos  b-cos  b'^  * 


Then  frcwi  (9)  we  see  that  uj^  (P)  is  given  by  (32)  with  C^^  replaced  by  C  , 
^ere 


/,r\    J,  _    i   (cos  c»cos  c')  sin  d 

K  "   Unk  (cos  a-coa  a'j(co3  b-cos  b'^  * 


Equation  (32)  for  any  u    is  identical  with  equation  (l5)  of  Part  I 
which  yields  the  field  on  rays  singly  diffracted  from  comers  of  the  edge.  Thus 
we  can  say  that  each  u    also  represents  the  field  on  such  rays,  and  that  C., 
C„  and  C-,  are  the  comer  diffraction  coefficients  ^ven  by  the  three  Kirchhoff 
theories.  These  coefficients  become  infinite  on  the  cones  a'  =»  a  and  b'  •  b. 
'i/hen  either  a'  ■  a  or  b'  ■  b,  an  edge  stationary  point  coincides  with  the  comer, 
and  a  different  method  of  evaluating  the  integrals  is  then  required,  llius  these 
cones  are  analogous  to  shadow  boundaries.  When  both  a*  -  a  and  b'  ■  bjan 
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Figure  U 
The  incident  ray  and  a  diffracted  ray  at  a  comer  of 
the  edge.  The  angles  between  the  incident  ray  and  the 
two  tangents  at  the  comer  are  a  and  bj  the  correspond- 
ing angles  for  the  diffracted  ray  are  a'  and  b». 
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•interior'  stationary  point  coincides  with  the  corner  and  a  still  different 

method  of  evaluation  is  required.  We  will  not  consider  either  of  these  cases 

further. 

Let  us' summarize  our  results.     We  have  seen  that  the  field  given  by  each  of 

the  three  ^irchhoff  methods  is,  for  large  k,   a  sum  of  three  parts.     One  part  is 

.lust  u         J  which  describes  the  incident  field  on  the  rays  which  pass  through 
"  geom 

the  aperture.     The  other  two  parts  describe,   respectively,  the  field  on  the  rays 

singly  diffracted  from  the  edge  and  on  the  rays  singly  diffracted  from  comejrs  of 

the  edge.     These  fields  have  exactly  the  form  determined  by  the  geometrical  theory 

of  diffraction.     Each  method  yields  an  edge  diffraction  coefficient  D  and  a 

comer  diffraction  coefficient  C  which  differ  from  the  coefficients  of  the 

geometrical  theory.     However,  the  diffracted  field  given  by  any  one  of  these 

methods  can  be  constructed  by  the  geometrical  theory  of  diffraction,  if  the 

appropriate  diffraction  coefficients  are  used. 

In  the  si^ecial  case  of  a  spherical  incident  wave,  u^  can  be  written  as 

u  plus  a  line  integral  along  the  aperture  edge.     This  was  shown  by 

geora  ^ 

Rubinowicz '    •* .     In  Appendix  III  we  show  that  the  asrrmptotic  evaluation  of  the 
line  integral  yields  exactly  u„     +  Uj,     ,  which  serves  as  a  check  on  our  previous 
calctilations. 


U.     Asymptotic  evaluation  of  the  Braunbek  integrals 

The  Bravinbek  approximations,  vriiich  we  have  called  Ug,   and  u-p*  are 
given  by  equations  (U)  and  (5)  respectively,  with  appropriate  integrands.     The 
integrands  are  determined  by  the  incident  field  at  points  away  from  the  edge,   and  by 
the  solutions  of  appropriate  half-plane  diffraction  probleiris  at  points  near  the 
edge.     We  will  assume  that  the  integrands  are  continuous  and  have  several  continuous 
derivatives.     Then  the  asymptotic  value  of  each  integral  for  k  large  comes  from 
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the  neighborhoods  of  the  three  types  of  critical  points.  Once  these  points  are 
found,  it  will  suffice  to  know  the  behavior  of  the  integrands  near  them  in 
oixier  to  evaluate  the  integrals  asjinptotically.  Therefore  it  will  not  be 
necessary  to  construct  a!q)licltly  the  complete  integrands,  which  join  the 
interior  and  edge  values  together  sraoothJ.y. 

In  the  interior  of  the  aperture  the  integrands  of  Ug,  and  Ugp  are  the 
same  as  those  of  u.  and  u^   respectively.  Therefore  u_,  and  u^jp  have  the  sane 
interior  stationary  points  as  u.  and  Up  and  also  receive  the  same  contribution 
from  each  point.  If  we  designate  the  sum  of  these  contributions  by  u^,  and  u^^ 
then  we  have  Ug,  "  ^  =  ^cecm*  ®"^  ^2  "  ^  "  ^seom*  ^^^  ^°*^  \  ^^'^  ^  equal 

^geom* 

II     II 
Let  us  now  locate  the  edge  stationary  points  and  evaluate  Ug,  and  Uggj 

the  contributions  to  Up,  and  yx^^   from  these  points.  To  do  this  we  must  consider 

the  integrands  near  the  edge.  We  will  express  them  in  terms  of  the  coordinates 

a-  and  s.  These  coordinates  are  unique  for  points  sufficiently  near  the  edge, 

provided  that  the  edge  has  no  comers.  Then,  according  to  Braunbek,  the  integrand 

at  (o-  ,s)  is  detemined  by  diffraction  from  a  half-plane  tangent  to  the  edge  at 

(o-,0).  The  field  incident  upon  the  half -plane  must  be  a  plane  wave.  Its 

amplitude,  phase  and  direction  must  coincide  with  those  of  the  actual  incident 

field  u.   ~  Ae^^^  at  the  point  (a-  ,0)  on  the  edge.  We  will  write  Mcr  ),   T(o-  ) 

and  VT(o-  )  to  indicate  the  values  of  A,  I  and  VT  at  this  point.  Then  the  plane 


wave  incident  upon  the  half-plane  may  be  written  as 

ik[T(o-)+r'VY(cr)] 


(37)     u^(r)  -  A(a')e 


Here  r  denotes  the  position  vector  of  a  point,  measured  from  an  origin  at  (o-  ,0) 
on  the  edge. 
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Now  let  u(cr  ,s)  and  u  (o" ,3)  denote  the  field  and  its  normal  derivative 

z 

at  the  point  (o-,s)  in  the  aperture  due  to  diffraction  of  the  plane  wave  (37) 
from  the  half -plane  tangent  to  the  edge  at  (cT  ,0),     Both  u  and  u     can  be 
detemdned  from  the  exact  solution  of  the  half-plane  problem^  ^  originally  solved 
by    A.  Sommerfeld.     In  using  this  solution,  it  is  necessary  to  take  account  of  tlie 
fact  that  VY(o-  )  makes  an  angle  p  with  the  edge  of  tlie  half-plane.     Then  the 
solution  for  the  screen  on  which  u  =  0  can  be  evaluated  in  the  aperture  and  yields 

(38)        u(cr,s)  -22±:^,(^)ei^t(-)-3n.VT(a-)l     (^^3) 

where 

fV.      a 


(39)        F^(o-,s)  = 


X2ks  sin  p  cos(|  -  |)         2 

e^^  dr. 


In  (38)  n  denotes  the  unit  normal  to  the  edge  pointing  into  the  aperture,  p 
is  the  angle  between  VT(<r  )  and  the  edge,  and  a  was  defined  previously  (see 
Figure  2),     Similarly  if  u-  ■  0  on  the  screen  we  obtain 

(UO)    u  (0-,  s)  -  2e:^  ,(,  )el^[Y(cr  )*   sn.VT(cr  )]  ^^   ) 

*^®^  y2ks  sin  p  C08(J  -  |)  _2 

(I4I)    F2(o-,3)  -  ik  sin  p  cos  o  e^'^  dV   -  sin(J  -  |)/^|p-^ 

2iks  sin  p  cos  (|-  -  ■j) 
We  now  insert  (38)  into  (U)  and  obtain  for  vu. 
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(U2)        ^l 


ke^"^     (    [  ,f   N^z    Uc[l(<r  )+sn.Vl(<r  )+rL  /         vi.     -1,    v    i ,   . 
-yn —  A(o-)-g-e      L^     '  JF^^Ccr  ,s)  |l-p     (<r)3|dsdo-. 

o     o 


Similarly,  inserting  (UO)  into  (^)  yields  for  tifg 

L 


ik[T(o'  )+sn«VT(a-  )+r] 


(U3)        u^  ~   "372~J     J  ^^^^^ S F,(or,s)|l-p'-^(o-)s|dsdo-. 


o     o 


In  (U2)  and  (U3)  the  upper  limit  of  the  s  integration  has  been  omitted  since, 
from  the  definitions  of  u_,  and  tigg*  only  contributions  fron  the  edge  at  the  lower 
limit  o-  ■  0  are  to  be  included. 

To  evaluate  (U2)  we  integrate  by  paHs  with  respect  to  a   and  obtain 


r  R 


ihh) 


Jo 


■"■  ikR[n.Vr+R  J 


-/. 


e^ 

HcR 


sn»VT*R 


h.W+R  ' 
s 


In  the  second  term  on  the  right  side  in  (UU)  we  have  retained  only  the  leading 

term  in  k,  which  comes  from  differentiating  F,  since  J7-   ■  0()^),  The  terms 

•»■      ds 

obtained  by  differentiating  other  factors  are  of  lower  order  in  k  since 

F.  ■  0(1),  and  therefore  they  are  omitted.  We  see  from  (39)  that  F-  «  0  at  s  "  0, 

so  the  first  tezm  on  the  right  side  of  (Ui)  is  zez*o.  Thus  when  __±  is  computed 

9s 

from  (39)  and  inserted,  (Ui)  becomes 
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(16)        /|2e^L»-^*\|l.p-^slds   ^    i/ipcos(5-|)   fe^'^'     '    '""    , 
//  ^  ^  ^       ^    Vo  >/i  R(R  +n.Vl) 


'  s 


where 


(li6)        I  -  R  +  3[n.VT+2  sin  p  cos^(J  -  |)] . 


To  evaluate  the  integral  in  {h$)  we  employ  the  following  asymptotic 


result'-     , 


(li7) 


gikh(s)     £(s2  dg    ^     A  g(0)  gikhCO+inA 


/S  ylch  (0~) 

s 


Before  using  this  in  (U5)  we  note  that  J,^  =  R  +n'VI+sin  p+sin  ^  sin  a.     At  s  ■■  0, 
R    "  VR.n  and  n«7T  =  -sin  p  sin  a.     Thus  5-  (<^)  "  sin  p+7R«n,     Now  using  (hi)  in 
ikS)  yields 

.  •    /I  /I       ci\o     ikR 

(W)       f  J|  «i>=[»n.mR]      |j     .13,^   „   le^  /;°ln6  '"^'t  '  t'V 

^  I  -jj  o  r^ii  H     a|ua  g        ^   2(.sin  p+VR.n;     (VT.n+VR.njR 


AH  quantities  in  the  right  member  of  (U8)  are  to  be  evaluated  at  s  -0,  i.e.  on 
the  edge.     Of  course  they  still  depend  upon  <r  . 

We  now  insert  the  right  side  of  (^8)  into  the  right-hand  member  of  (U2), 
which  becomes  a  single  integral  with  respect  to  <r  ,     This  integral  can  be 
evpluated  by  stationary  phase.     In  fact  the  phase  of  the  integrand  is  just  R+T, 
which  is  exactly  the  same  as  the  phase  of  the  integral  for  a.      in  equstion  (20), 
Therefore  the  stationary  points  are  the  same  as  those  previously  considered  and 
the  result  of  the  evaluation  is  again  of  the  form  (21)  with  TL    replaced  by 
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(U9) 


.^"\c03(f  .  \) 


Bl 


ynk  sin  pCsin  p+VR«n)(VYT^VR.n) 


In  terms  of  the  angles  previously  defined  we  have  R  ■  sin  p  cos  9, 
Vr»n  =  -sin  p  sin  o,  and  VR»n  =  -sin  p  sin  0,  Upon  substituting  these  ejrpressions 
into  {hS)   and  simplifying  the  result,  we  obtain 


(5b) 


e 


inA 


^Bl 


2/2nk  sin  p 


sec  -xC©  -  a)+csc  --(e  +  o) 


Now 


II 


tL^  is  given  by  (21)  with  Dg,  inserted  in  place  of  D.  .  The  resulting 


expression  is  eocactly  the  same  as  the  field  given  by  the  geometrical  theory  of 
diffraction  on  rays  singly  diffracted  from  edges.  Therefore  u^^  can  also  be 
interpreted  as  the  field  on  the  edge  diffracted  rays  and  IL^  as  the  edge 
diffraction  coefficient  given  by  the  Braunbek  method  for  screens  on  which  u  »»  0. 
This  edge  diffraction  coefficient  is  exactly  the  same  as  that  of  the  geometrical 
theory  for  screens  on  which  u  =  0,  which  is  given  in  Part  I,  equation  (13 )• 


n 


A  procedure  similar  to  that  used  for  u^,  applies  to  the  evaluation 


II 


of  u^o.  We  first  insert  (Uo)  into  (5)  and  obtain 

L 


(51)    4^ 


'^ri/h   r  f             ^ik[T(o-)+sn.Vl(or  )+r]        ,    , 
"2-^^    A({r)2 iF2(cr,3)|l-p'-'(a-)s|dsd<r 


o  o 


Next  we  integrate  by  parts  with  respect  to  s,  but  we  only  retain  the  fii^t  term 
in  F_(cr,s)  since  that  is  the  leading  teiw  in  k.  In  this  way  we  get  an  integral 
of  the  type  given  in  (U?).  Then  using  (U?)  we  finally  obtain 
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(52) 


f 


A((r 


.ik  I+sn.VT+R 


e 


It    ^-1    I J  iTi/lii/n   gin  6 

|l-p  sids  ~  -e   /  ^(3i„  LyR.nj 


n   a> 


sin  p  cos  a  cos(i —  •=) 


Vf«n  +  7R»n 


+  sin(^--j) 


ikR 


II 


Now  inserting  (52)  into  (U3)>  we  express  Ur.r,   as  a  single  integral  with  respect 
to  o"  ,  Evaluating  this  by  stationary  phase  again  leads  to  (21)  with  D,  replaced 


by  Dg2  which  is 


.inA 


-   ^.=^;Q^ 


cos  a  cos('jj  -  Tj)    s^(r  -  "j) 


VJ«n  +  VR»ii 


sin  p 


In  terms  of  the  angles  previously  considered  this  becomes 


(5U) 


D 


II 
B2 


.in  A 


2;^iarsin  p 


sec  •»(©-  a)  -  CSC  "k(©+  a) 


This  Braunbek  diffraction  coefficient  is  exactly  the  same  as  that  given  by 
the  geOTietrical  theory  of  diffraction,  Part  I, equation  (13) »  for  a  screen  on 
which  the  normal  derivative  of  the  field  is  zero.  Thus  for  such  screens  the 
Bravinbek  method  also  gives  exactly  the  same  field  on  rays  singly  diffracted 
from  the  edge  as  does  the  geometrical  theory. 

If  we  attempt  to  determine  u_^  or  Up«  ,  the  contributions  from   the 
comers  of  the  edge,  we  find  that  the  Braunbek  method  fails.  This  is  due  to  the 
fact  tliat  a  point  in  the  aperture  near  a  comer  is  near  two  different  points  of 
the  edge.  Therefore  the  Braunbek  method  yields  two  different  values  for  the 
integrand  at  such  a  point.  Consequently  the  integrand  near  a  comer  is  not 
uniquely  determined,  so  the  comer  contribution  cannot  be  evaluated.  However, 
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the  Braunbek  method  can  be  modified  by  assuming  that  the  field  near  a  comer 
is  the  same  as  that  near  an  infinite  comer  hit  by  a  plane  wave.  By  an  infinite 
comer  we  mean  a  screen  bovinded  by  two  semi-infinite  straight  lines  which  meet 
at  the  comer.  This  modification  irould  undoubtedly  yield  a  field  on  comer 
diffracted  rays  in  agreement  with  that  given  by  the  geometrical  theory  of 
diffraction.  But  this  cannot  be  verified  at  present  because  of  the  intractable 
form  of  the  solution  of  the  comer  diffraction  problem. 

In  summary,  we  have  found  that  the  Braunbek  method  applies  to  an  aperture 
without  comers.  For  large  k  it  yields  a  field  which  is  a  sum  of  two  parts. 
One  is  ii^---*  the  incident  field  on  the  rays  which  pass  through  the  aperture. 
The  other  is  the  field  on  rays  singly  diffracted  from  the  edge.  This  field  is 
of  exactly  the  same  form  as  that  given  by  the  geometirLcal  theory  of  diffraction. 
For  the  two  types  of  screen  considered  the  diffraction  coefficients  are  precisely 
the  same  as  those  given  by  the  geometrical  theory.  These  coefficients  are 
compared  with  the  Kirchhoff  coefficients  in  Figure  ^  .  However  the  Braunbek  method 
does  not  take  account  of  multiple  diffraction. 

At  caustics  and  on  shadow  boundaries  the  evaluations  of  the  Braunbek 
integrals  which  we  have  given  are  invalid.  At  such  places  the  integrals  can  be 
evaluated  by  a  slight  modification  of  the  present  procedure, 

5.  The  transmission  cross  section 

The  transmission  cross  section  cr  of  an  aperture  is  defined  as  the  ratio 
of  the  flux  through  the  aperture  to  the  flux  per  unit  area  in  the  incident  wave, 
assuming  that  the  incident  wave  is  plane.  According  to  the  cross  section 
theorem,  o-  can  be  determined  from  the  amplitude  of  the  transmitted  field  in 
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the  forward  direction,  i.e.  in  the  direction  of  propagation  of  the  incident 
wave.  Let  the  field  in  this  direction  at  a  large  distance  R  from  the  aperture 
be  given  by 

,  ikR 

Here  f  is  a  constant.  Then  if  the  phase  of  the  incident  wave  is  zero  at  the 
origin  of  R,  the  cross  section  theorem  states'-  -* 

(56)  o-  -  Im  f  . 

Of  course,  this  theorem  applies  to  the  exact  solution  u.  However  we  will  apply 
it  to  the  approximate  solutions  which  we  have  been  considering. 

Let  us  first  assume  that  the  incident  wave  is  e^^^  sin  y  +  z  cos  y)  ^ 
Then  (7)  jrields 

T  r   ^ik(x  sin  f  +  z   cos  y  +  R)       1 

(57)  u,(P)  -  ^^j^  2 R^  (Ik  -  l)d»J7  . 


If  P  is  far  from  the  aperture  in  the  forward  direction,  the  phase  in  (5?)  is 

essentially  constant  over  the  aperture.  Furthermore  R  is  practically  -cos  k 

z 

and  1/r  is  negligible  compared  to  ik.  Thus  (5?)  yields 


ikR 
(58)     ti^CP)  ~  -  ^^   .  lA  cos  Y 


Here  A  is  the  area  of  the  hole.  Now  (56)  yields 
(59)     0-,  =  A  cos  Y  . 
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Tife  see  that  the  first  modified  Kirchhoff  approximation  yields  for  O"  exactly 
the  projected  area  of  the  aperture.  This  is  just  the  geonetrical  optica  cross 
section.  This  result  is  precise  and  not  merely  valid  for  large  k.  The  second 
modified  Kirchhoff  approximation  jrLelds  the  same  result  and  therefore  so  does 
the  original  Kirchhoff  method. 

If  we  now  consider  the  Braunbek  method,  we  find  that  the  interior  points 
of  the  aperture  also  contribute  the  expression  in  (58)  to  the  field  and  A  cos  y 
to  the  cross  section.  This  is  so  because  at  interior  points  the  Braunbek  integrand 
is  the  same  as  the  Kirchhoff  integrands  wiiich  we  have  Just  considered.  However 
there  is  an  additional  contribution  to  the  Braunbek  field  from  the  neighborhood 
of  the  edge.  This  is  given  in  equations  (U2)  and  (U3)  with  A(o- )  »  1,  ¥  =  x  sin  y 
and  Vf  =  (sin  y.  0»  cos  y)»  When  these  values  are  inserted  the  phase  of  the 
integrand  becomes  constant  and  equal  to  R  at  points  far  from  the  aperture  in  the 
forward  direction.  Then  R  is  practically  equal  to  -cos  y  ^^   R  is  essentially 
constant  over  the  aperture.  Then  (li2)  and  (U3)  become 

(60)     ug  -"^^^3/2"  °°'^  [  f  F^(<r,s)|l-p-^sldsdcr  , 


-.-.     ikR-in/li  ^     r  , 

(61)  ujg  -  "^-775 F2(o-,s)ll-p"-^sId3dcr  . 

'  o  o 

To  evaluate  the  integral  in  (60)  we  first  integrate  by  parts  with 
respect  to  s,  ignoring  the  contribution  from  the  upper  limit,  and  obtain 

(62)  rF^(o  ,s)|l-p'^s|d3  ^  -  J  )^({2s^^^   -  p"y/^)e^^ds  . 
^  n  y  a 


-2$- 


Here  we  have  introduced  p,  defined  by 


(63)     p  »  2  sin  p  cos^(J  -  |) 


Now  integrating  the  first  term  on  the  right  of  (62)  by  parts  once,  and  the 
second  teim  twice,  we  get 


(6Ii) 


|F^(<r,s)|l-p-^|ds    ^.M. 


i    +       3 

,  Uc  p  P 


X 


^ikps 


/i 


ds 


u 


^  i;?7pi 


n  ^inA 

T5^ 


In  the  second  line  of  (6U)  we  have  inserted  the  asymptotic  value  of 


I 


3-1/2  ^iicps  ^^ 


Upon  inserting  (6U)  into  (60),  and  adding  to  this  (58),  which  represents  the 
contribution  of  interior  points,  we  obtain 


(65)   ^(p)„.!gf[ucosr*^J  %■ 


-  3i  cos  Y  [  do;: 
Sk^    J^  P  P 


Using  the  cross  section  theorem  we  now  find 


(66)    cr  ^,  ^  A  cos  Y  -  ^-2^ 


KL 


8k^ 


<io- 


In  the  case  of  normal  incidence  (y  =  0)  (65)  and  (66)  become 
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(67)    V^)  -    ^[^*h'^]> 


(68)   o-g^   ^  A  -i2^  . 

Uk 

The  evaluation  of  iLp  is  practically  the  same  as  that  of  vl, .  The 
integrand  now  involves  two  terms  since  F^  is  a  sum  of  two  tenns.  The  first  of 
these  terms,  which  is  the  only  one  we  retained  in  our  previous  calculations, 
is  just  ik  sin  a  cos  p  F-(a,s).  Since  a  and  p  depend  upon  cr  only,  this  term 
yields  the  second  two  terms  in  (65),  but  with  the  factor  cos  6 
replaced  by  sin  a  cos  p  under  the  integral.  Therefore  in  place  of  (65)  and  (66) 
we  obtain 

L  L 


ikR 

(69)   u^W  ~  =^ 


r 


iA  cos  6  •.  ^    «^  P  "°^  °  d^  -  ^    ^^gP  "°"  °  do- 
K  8k  j^   p  p 


(70)    o-„„  ^  A  cos  6  -  -2^    ^^  E  ^°^  ^  do-  . 


At  normal  incidence  (69)  and  (70)  yield  (6?)  and  (68)  for  ^^^(P)  a«i  ''"b2*  ^®° 

the  second  term  in  the  integrand  of  u^p  is  evaluated,  it  adds  terms  in  l/1c  and 

2  2 

i^  to  the  bracketed  expression  in  (69)  and  a  tenn  in  lA  "to  (70), 

Let  us  now  coitpare  the  results  of  this  section  for  the  transmission  cross 
section  with  other  results.  Consider  first  the  case  of  normal  incidence  on  a 

circtJ.ar  hole  in  a  screen  on  which  u  »  0 .  In  Part  I  it  is  shown  that  the  first 

-3/2 
correction  to  the  geometrical  cross  section  in  this  case  is  of  the  order  k  '  . 

This  correction  is  obtained  by  taking  account  of  doubly  diffracted  rays.  The 

Kirchhoff  methods  yield  no  correction  and  the  Brauribek  method  yields  (68),  all 
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failing  to  yield  the  k  correction.     Next  consider  arbitrary  incidence  upon 

an  aperture  of  any  shape  in  a  screen  on  which  u  =  0,     In  Part  I  it  is  shown  that 
usually  the  doubly  diffracted  rays  lead  to  a  correction  of  order  k      in  this 
case,  just  as  the  Brauribek  result  (68)  does.     However  the  result  based  upon 
doubly  diffracted  rays  involves  the  aperture  size  and  shape,  while  the  Braunbek 
method  yields  the  same  correction  for  all  apertures,     Finally  let  us  consider 
noiwal  incidence  on  a  circular  hole  in  a  screen  on  which  the  normal  derivative 
of  u  is  zero.     In  this  case  H,  Levine^  ■'   has  shown  that  the  correction  is  -  *— w 

instead  of  -  ^-k,  which  is  given  by  the  Braunbek  method  in  (68),     VJhen  the 

^  II 

second  term  in  the  Braunbek  integrand  for  ulp  is  included  the  disagreement  ia 

further  increased. 

On  the  basis  of  these  comparisons  we  must  conclude  that  the  Braunbek 

method  fails  to  yield  the  proper  correction  to  the  cross  section  when  the  cross 

section  theorem  is  used.  However  the  Braunbek  method  does  yield  the  correct 

field  on  singly  diffracted  rays,  since  it  agrees  with  the  field  given  by  the 

geometrical  theoiy  of  diffraction.  The  disagreement  in  the  former  case  and 

the  agreement  in  the  latter  case  can  be  explained  as  follows.  In  calciilating 

the  field  on  singly  diffracted  rays  use  was  made  of  only  the  first  two  derivatives  of 

the  Sommerfeld  half-plane  solution  along  the  nonnal  to  the  edge,  evaluated  at  the 

edge.  On  the  other  hand,  in  computing  the  cross  section  correction  additional 

derivatives  of  the  Sommerfeld  half-plane  solution  were  used.  We  may  therefore 

conclude  that  only  the  first  two  derivatives  of  the  aperture  field  at  the  edge 

are  correctly  given  by  the  half-plane  solution  for  k  large.  This,  then,  is 

perhaps  the  proper  statement  of  the  Braunbek  hypothesis. 


6»  Conclusion 

We  have  evaluated  asymptotically,  for  small  \  or  large  k,  the  integrals 
which  give  the  field  diffracted  by  an  aperture  of  any  shape  in  a  thin  screen. 
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The  five  integrals  which  we  have  considered  are  those  given  by  the  Kirchhoff 
method,  its  two  usual  modifications,  and  Braunbek's  new  modification  for  two 
types  of  screen.  In  every  case  the  asymptotic  form  of  the  field  is  a  sum  of 
contributions  from  three  types  of  critical  points.  These  are  interior  stationary 
points,  edge  stationary  points,  and  comers  of  the  edge.  All  the  methods  lead 
to  exactly  the  same  critical  points,  and  these  points  correspond  precisely  to 
certain  rays  given  by  the  geometrical  theory  of  diffraction.  Interior  stationary 
points  correspond  to  rays  directly  transmitted  through  the  aperturej  edge 
stationary  points  correspond  to  rays  singly  diffracted  at  the  edge,  and  comers 
correspond  to  rays  singly  diffracted  from  comers. 

The  contribution  to  the  field  from  each  type  of  point  is,  for  every 
method,  of  exactly  the  form  given  by  the  gecsnetrical  theory  of  diffraction  for 
the  corresponding  type  of  ray.  Each  method  yields  certain  edge  and  comer 
diffraction  coefficients,  except  the  Brauribek  method,  idiich  does  not  apply  to 
apertures  with  comers.  The  edge  diffraction  coefficients  given  by  the  Brauhbek 
method  are  preciselj''  the  same  as  those  given  by  the  geometrical  theory  of 
diffraction  for  both  types  of  screen,  Ihus  the  Braunbek  method  yields  exactly 
the  same  field  as  that  given  by  the  direct  and  singly  edge  diffracted  rays  of 
the  geometrical  theory  of  diffraction.  It  could  be  modified  to  apply  to 

comers  but  it  takes  no  account  of  multiple  diffraction,  as  does  the 
geometrical  theory.  The  other  Kirchhoff  methods  yield  fields  which 

correspond  to  direct  rays,  to  singly  edge -diffracted  rays  and  to  singly  comer- 
diffracted  rays.  They  yield  the  same  field  on  the  direct  rays  as  do  the 
geometrical  theory  and  the  Bravinbek  method.  However  they  yield  fields  on  the 
edge  diffracted  rays  which  differ  from  each  other  and  from  the  Braunbek  method, 
which  agrees  with  the  geometrical  theory. 

We  conclude  that  the  resvilts  just  described  lend  considerable  support  to 
the  geometrical  theory  of  diffraction.  They  also  show  that,  for  apertures  without 
comers,  the  Braunbek  method  is  equivalent  to  the  geometrical  theory  provided 
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that  only  singly  diffracted  rays  are  considered.  Since  the  geometrical  theory 
is  simple  to  use  and  requires  no  integration,  it  seems  to  be  preferable  to  the 
Braunbek  method.  It  can  also  take  account  of  comers  and  of  multiple  diffraction. 
On  the  other  hand,  the  geometrical  theory  involves  supplementary  considerations  in 
oixier  to  describe  the  field  on  caustics  and  shadow  boundaries.  The  Braunbek 
method  encounters  no  difficulty  on  such  surfaces.  Therefore  it  can  be  used  to 
supplement  the  geometrical  theoiy. 

We  also  see  why  the  simpler  Kirchhoff  methods  predict  the  correct 
qualitative  behavior  of  the  diffracted  field.  This  is  because  the  Kirchhoff 
methods  yield  the  sacie  singly  diffracted  rays  and  the  same  phase  as  do  the  more 
accurate  geometrical  theory  and  the  Braunbek  method.  Since  the  qualitative 
properties  of  the  diffraction  pattern  depend  upon  interference  between  the  fields 
on  these  rays,  they  are  predicted  essentially  correctly.  Quantitatively  the 
Kirchhoff  methods  are  also  fairly  accurate  in  and  near  the  forward  direction. 
This  is  because  the  Kirchhoff  diffraction  coefficients  are  close  to  the  more 
correct  coefficients  of  the  Braiinbek  method  and  of  the  geometrical  theory  of 
diffraction  in  and  near  the  forward  direction.  However  away  from  this  direction 
the  Kirchhoff  methods  must  be  quantitatively  inaccurate.  They  will  also  fail 
at  smaller  values  of  k  when  multiple  diffraction  becomes  important. 

Finally  we  have  seen  that  the  Kirchhoff  method  yields  no  correction^  and 
the  Braunbeck  method  yields  an  improper  correction,  to  the  geometrical  optics 
transmission  cross  section  of  an  aperture  when  it  is  computed  by  means  of  the 
cross  section  theorem.  For  a  screen  on  which  u  "  0,  this  is  due  to  ths  failure 
of  these  methods  to  take  account  of  multiple  diffraction.  In  the  geometrical 
theory  of  diffraction  the  appropriate  correction  to  the  cross  section  for  a 
circular  hole  in  such  a  screen  is  obtained  as  soon  as  doubly  diffracted  rays 
are  taken  into  account* 
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Appendix  I.     Interpretation  of  u 

We  have  aastmed  that  the  incident  field  is  asymptotically  (for  large  k) 

ikT  Fl  nl 

of  the  form  Ae       .     It  then  follov/s  •■    -^  from  the  wave  equation  that  T  satisfies 

the  eiconal  equation  (VY)     -  1  and  that  A  is  given  by 
(Al)  A(P)   = 


y(l+p~^R)(l+P2-4l) 


In  (Al)  A  denotes  the  value  of  A  at  some  point  P  on  a  ray,  A(P)  is  the  value 

of  A  at  P,  a  distance  R  along  the  ray  fran  P^,  and  p^  and  p^   are  the  principal 

radii  of  curvature  of  the  wavefront  T  =  const. at  P  •  The  result  (Al)  also 

follows  from  the  law  of  conservation  of  energy  within  a  tube  of  rays.  From  the 

eiconal  equation  we  find  that  T(P)  =  ^q+R* 

We  wish  to  show  that  each  terra  in  equation  (1^)  for  u^(P)  and  Ti2(P)  is 

of  the  form  A(P)e      where  A(P)  and  T(P)  are  deteimined  as  indicated  above. 

It  will  then  follow  that  both  u^  (P)  and  U2(P)  equal  ^^^0^^^*  if  P  is  in  the  lit 

region.  From  (l^)  we  see  at  once  that  the  phase  has  the  desired  form  and  that 

the  amplitude  is  proportional  to  the  value  of  A  at  the  stationary  point.  Therefore 

calling  this  point  P  ,  xre  may  identify  A(P  )  with  A  in  equation  (Al)  above.  It 

remains  to  be  shown  that  the  coefficient  of  A(P^)  in  (15)  equals  that  of  A^  in 

equation  (Al), 

The  factors  in  the  denominator  of  (Al)  may  be  negative  if  the  radii 

p^  and/or  p^   are  negative.  If  one  factor  is  negative  the  root  to  be  taken  is 

e  '  multiplied  by  a  positive  real  number.  If  both  factors  are  negative  the 

root  to  be  talien  is  e   multiplied  by  a  positive  real  number.  These  choices 

correspond  to  the  fact  that  the  phase  on  a  ray  changes  by  -  ■«  each  time  the 

ray  goes  through  a  caustic'-  •'  and  when  this  occurs  a  factor  of  the  denominator 

becomes  negative.  Therefore  we  can  rewrite  (Al)  by  introducing  e'  which  is  +1, 

-l,or  -i  according  as  neither,  both, or  one  center  of  cuivature  of  the  wavefront 

lies  between  P^  and  P.  If  we  also  introduce  the  mean  curvature  M  ■=  fC  +  p" 
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and  the  Gaussian  cuivature  C  =  Pi'Pj'''  o^  ^^^  wavefront,  (Al)  becomes 

e'A^ 

(A2)    A(P)  =.  -  ■■   . 

/  ll+MR+CR'^1 

This  must  be  compared  with  the  coefficient  of  A(P^)  in  (1^),  which  is 
——————  p     p 

e  /(1-y2.I^)/|AI.  To  evaluate  A  "«  note  that  R^  -  f^:^     =  ^i;?2^=  i^  , 


id-th  similar  expressions  for  IL.^  and  R  .  Inserting  these  eoqsressions  into  the 
definition  oH  /\,   equation  (l6),  the  coefficient  of  A(P^)  becomes 


>    i-il-il  1  -  "^  -  V  - 


It  is  now  easy  to  show  that  the  coefficient  of  IT  in  (A3)  is  just 

the  Gaussian  curvature  C  and  the  coefficient  of  R  is  the  mean  curvature  M  of 

the  wavefi^nt  at  P  .  It  also  follows  from  their  definitions  that  e  and  e' 
o 

ai«  equal.  Thus  the  coefficient  of  A  in  equation  (A2)  is  equal  to  the 
coefficient  of  A(P  )  in  (l5)  which  is  given  by  equation  (A3).  This  completes 
the  proof  of  the  fact  that  both  u^(P)  ani  ij^(P)  are  equal  to  ^ijjj.(P)*  if  P  is 
in  the  lit  region. 

Appendix  II.  The  geometrical  factor  G 

In  order  to  siniiLify  the  expression  for  G  given  in  (22)  we  must 
evaluate  (Y+R)^  ^  at  an  edge  stationary  point.  To  do  so  it  is  convenient  to 
introduce  t,  the  xrnit  vector  tangent  to  the  edge.  Then  (T+R)   =  VT»t  +  7R»t  " 
cos  p  +  TR't,  Therefore  (T+R)^.  ^  »  -^^  sin  p  +  VR»t  +  (VR)  't  ■  -p  sin  p 

♦p'^^'n  +  R"  sin*^.  Here  we  have  substituted  t   =  p  n  where  p  is  the  radius  of 

—1   2 

curvatura  of  the  edge  and  n  the  unit  normal.  In  showing  that  (VR)  't  »  R  sin  p 

WB  made  use  of  the  fact  that  7R«t  *  -cos  p  at  the  stationary  point  and  we  perfonned 
a  sin^^le  calculation.  If  we  now  define  the  angle  5  by  cos  6  ■  -VR»n  then  (22) 
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becomes 


(All)         G  -  g^P  e"^"^  *  ^"^ 

R/l-P^y.sin  p-pcos&fRsin^P  I 

Let  us  now  define  p,  by 

(A5)        p,  = =5i5^^:3 . 

p^sin  p  +  p    cos  6 


Then  (Ah)  becomes 

(A6)    G  =  j|R(l+p^^|]'^^^  e-^"^  *  ^"^ 


If  we  omit  the  exponential  factor  and  the  absolute  value  signs  under  the  radical, 
G  is  unchanged  provided  that  the  square  root  is  determined  in  the  manner  explained 
just  preceding  equation  (A2),  This  is  obviously  the  case  if  the  radicand  is 

—in  /? 

positive,  while  If  it  is  negative  the  root  introduces  the  necessary  e"  '  . 
Thus  we  obtain  equation  (26),  as  we  desired  to  show.  It  is  to  be  noted  that 
the  definition  of  p.  above  in  (a5)  is  exactly  the  same  as  that  in  equation  (11) 
of  Part  I. 
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AppendJJC  III.  Asymptotic  evaluation  of  Rublnowicz' line  integral 

Suppose  that  the  incident  wave  is  the  spherical  wave  e   /r,  where  r 

denotes  distance  from  the  source.  A,  Rubinowicz'-  J  has  shown  that  then  u„  can 
be  written  in  the  fom 

(n)        „  -„      1  r  e^^'^^^  sin  p  m«VR  ., 


\  '    geom  '  Hn  j   iriR      l+Vr.VR 


The  integral  is  a  line  integral  along  the  edge  of  the  aperture  and  the  vector  m 
is  the  unit  outward  normal  to  the  cone  formed  by  lines  joining  the  source  to  the 
points  of  the  edge.  Rubinowicz  has  evaluated  this  integral  by  the  method  of 
stationary  phase.  We  find  that  this  procedure  applies  to  edge  stationary  points 
and  yields 


(A8)    uj^  ~  Z  Dj^r-^e^^^^^G 


Here  the  sum  extends  over  the  edge  stationary  points,  >rtiich  are  determined  by 
the  condition  (r+R)  «  0.  These  are  Just  the  same  stationary  points  as  those 
obtained  in  evaluating  (20),  The  factor  G  is  given  by  equation  (22)  and  D. 
is  defined  by 


U9)       v,..^^±S^ 

2  >^iS  (l+Vr.VR) 


From  its  definition,  we  see  that  m  is  perpendicxilar  to  both  Vr  and  t. 
Since  it  is  a  unit  vector,  we  have  m  «  tJCVp/sin  p.     Then  m«VR  ■  t3cVr*VR/sin  p. 
In  order  to  evaluate  this  triple  product  it  is  convenient  to  consider  the  normal 
vector  n,  pointing  into  the  aperture,  and  the  binormal  b,  pointing  toward  the 
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dai4c  side.  Then  we  have 

(AlO)   Vr  »  t  cos  p  -  n  sin  p  sin  a  +  b  sin  p  coa  a  , 

(All)   VR  -  -t  cos  p  -  n  sin  p  sin  9  +  b  sin  p  cos  9  . 

Upon  evaluating  the  triple  product  and  the  scalar  product  Vr«VR,  and  inserting 
the  results  into  (A9)  we  obtain 

/I 
(A12)      D  ■   e  '^  (sin  g  coa  9  -  cos  a  ain  9) 

2  yS±   sin  p(l+sin  a  sin  9+cos  a  cos  9) 

After  some  trigonometric  transformations  we  find  that  IL  given  by  (A12)  is 
exactly  the  same  as  Djj  given  by  (29).  Therefore  the  stationary  points  of  the 
line  integral  contribute  exactly  the  game  terms  to  the  field  vl.   as  do  the  edge 
■tationary  points  of  the  double  integral  (6),  which  also  represents  tl.. 

Let  us  now  consider  the  contributions  of  comers  to  the  line  integral. 
These  contributions  can  be  obtained  by  integrating  by  partsi  and  ignoring  the 
2>emaining  integral  which  is  of  lower  order  in  k.  In  this  way  we  obtain 


III         1   eil«(r^R) 
(A13)   u^   ~  Z  -  HT"  rfeiJc(l^l>r.rRJ 


sin  p  m»VR 


o- 


The  brackets  indicate  the  jtunp  or  discontinuity  of  the  enclosed  quantity  at  a 
comer,  and  the  sum  extends  over  all  comers.     This  sum  can  be  rewritten  in  the 
forai 

TTT  1       .ik(r+R) 

(All.)       u^^    ^      rc,.i.5_X—    > 
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where 


/.-,r'^    r.        i      r  sin  p  m«VR  1 


Upon  inserting  the  angles  previously  defined,  we  find  that  (Al5)  becomes 
exactly  (36). 
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